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The origin of x-ray magnetic circular dichroism (XMCD) at the Cu 1/2,3 edge in 
YBa2Cu307/Lai_ a; Ca I Mn03 superlattices is revealed by performing first-principle electronic struc- 
ture calculation using fully-relativistic spin-polarized linear muffin-tin orbital and projected aug- 
mented plane wave methods. We show that the XMCD spectra at the Cu L2,3 edges are proportional 
to the difference of the densities of Cu d^ 3z 2_ 1 and d^ 3z 2_ 1 states. Although the Cu d 3z 2_ 1 states lie 
well below the Fermi level, a small number of d^ 3z 2_ 1 holes is created by the Cu d 3z 2_ 1 -OBa p z -Mn 
d 3z 2_ 1 hybridization across the interface. Even this tiny number of holes is sufficient to produce 
appreciable Cu L2,3 XMCD. The robustness of this conclusion is verified by studying the influence 
of doping, atomic relaxation, correlation effects, and antiferromagnetic order in a Cu02 plane on 
the XMCD spectra. 



PACS numbers: 75.50.Cc, 71.20.Lp, 71.15.Rf 
I. INTRODUCTION 

Artificially grown heterostructures and superlattices 
of alternating superconducting (SC) and ferromagnetic 
(FM) materials have become an important tool for ex- 
ploring the interplay between superconductivity and fer- 
romagnetism. Various interaction phenomena that may 
occur, like proximity effects, coupling phenomena, charge 
transfer etc., are directly influenced by the structure and 
especially by the physical and/or chemical disorder of the 
interfaces that separate them. The combination of ferro- 
magnetic and superconducting materials could end up in 
devices like a spin-controlled transistor with a high gain 
current and short switching times. 

In the last years much attention has been payed 
to oxide based FM/SC superlattices (SL), combin- 
ing a high-T c superconductor YBa2Cu3C>7 (YBCO) 
with a La2/3Cai/3Mn03 (LCMO) manganite which ex- 
hibits colossal magnetoresistancej 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 Both 
YBCO and LCMO have oxide perovskite structure with 
very similar in-plane lattice parameters, which allows 
the growth of superlattices with sharp interfaces, thus 
strongly reducing extrinsic " structural" effects which oth- 
erwise could obscure the FM/SC interplay. Using scan- 
ning transmission electron microscopy and high spa- 
tial resolution electron energy loss spectroscopy, Varela 
et al. showed in Ref. @ that individual layers of the 
YBa 2 Cu307/Lai_ a; Ca 3 ;Mn03 superlattices are flat over 
long lateral distances. The interfaces are coherent, free of 
defects, exhibiting no roughness, and are located at the 
BaO plane of the superconductor. Concerning chemical 
disorder, EELS measurements show the absence of mea- 
surable chemical intcrdiffusion within experimental error 
bars. 

The work on the perovskite oxide FM/SC superlattices 
is motivated by the appealing properties of the cuprate 
high-T c superconductors (HTSC) whose high SC critical 



temperatures make them potentially useful for technolog- 
ical applications. Further, since HTSC are believed to be 
susceptible to a variety of competing instabilities, there 
is a high potential for SC/FM quantum states in multi- 
layer structures. It has been suggested by Si in Ref. [H 
that electron injection from hole-doped rare earth man- 
ganites into high-temperature superconductors can yield 
information on spin-charge separation in high-T c super- 
conductors. The spin-polarization, in addition to the 
ferromagnetic ordering nature, is one particular prop- 
erty of hole-doped rare earth i2i_ x j4 K Mn03 manganite, 
where R is a trivalent rare-earth ion and A is a diva- 
lent alkaline-earth ion. In the case of La 2 /3Ca! /3M11O3 
a nearly full spin polarization of the transport electrons 
can be founds 

Magnetization measurements of epitaxially grown bi- 
layers of LCMO and YBCO show a coexistence of ferro- 
magnetism and superconductivity in these bilayer sam- 
ples at low temperatures. 1 - Both magnetization measure- 
ments which are used to determine T c and magneto- 
optical measurements for the j c evaluation of the su- 
perconducting film show up a critical thickness for the 
YBCO film of about d=30 nm, below which the su- 
perconductivity is strongly affected by the magnetic 
film. Sefrioui et al. in Ref. 3 studied the interplay be- 
tween magnetism and superconductivity in high-quality 
YBCO/LCMO superlattices. They found evidences for 
the YBCO superconductivity depression in the presence 
of the LCMO layers. Due to its short coherence length, 
superconductivity survives in the YBCO down to a much 
smaller thickness in the presence of the magnetic layer 
than in low T c superconductors. For a fixed thickness 
of the superconducting layer, superconductivity is de- 
pressed over a thickness interval of the magnetic layer 
in the 100 nm range. 

Holdcn et al. in Ref. reported ellipsometry measure- 
ments of the far-infrared dielectric properties of SL com- 
posed of thin layers of YBCO and LCMO. The optical 
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data provide clear evidence that the free-carrier response 
is strongly suppressed in these SL as compared to the 
one in the pure YBCO and LCMO films. The suppres- 
sion occurs in the normal as well as in the SC state and 
involves a quite large length scale. 

This research is in its early stage, and relatively little 
is known about the nature of magnetism at the interface, 
the spatial distribution of the magnetization throughout 
the layers, and the interplay of FM and SC order param- 
eters in general. 

X-ray magnetic circular dichroism (XMCD) spectra at 
Mn and Cu L 2 ,s edges in YBCO/LCMO SL were mea- 
sured by Chakhalian et al. in Ref. ITol . It was shown that 
magnetic dichroism is clearly present at both Mn and 
Cu edges, although the Cu XMCD signal is small com- 
pared to that of Mn (27%) and does not exceed 1.4%. 
The measurements suggest the presence of an uncompen- 
sated induced magnetic moment in a YBCO layer close 
to the YBCO/LCMO interface. As manganite layers un- 
dergo a ferromagnetic transition at around 180 K, the 
large dichroism at the Mn £2,3 edge is expected. How- 
ever, the presence of the net ferromagnetic magnetiza- 
tion on Cu is quite surprising. Authors investigated also 
the temperature dependence of the dichroic signals at 
both the Mn and Cu £2,3 edges. At temperatures above 
170 K, the dichroism on Mn disappears, consistent with 
the ferromagnet-paramagnet transition. Similar behav- 
ior was also found for Cu, except that already above 150 
K its XMCD signal decreases below the detection limit. 
From the similar temperature dependences of the XMCD 
spectra authors concluded that the magnetic moment on 
Cu is induced by strong interactions between Cu spins 
and the ferromagnetic moment of Mn across the inter- 
face. 

Later, x-ray absorption spectra (XAS) and x-ray linear 
dichroism (XLD) at Mn and Cu £2,3 edges were recorded 
in the total electron yield (TEY) and fluorescence-yield 
(FY) modes from the same sample of YBCO/LCMO 
SLfii In contrast to strong Cu £2,3 XLD detected in the 
bulk sensitive FY mode, only a very weak XLD signal was 
measured in the TEY mode, which is supposed to probe 
interfacial Cu ions. This led authors to the conclusion 
that the electronic states of Cu ions at the YBCO/LCMO 
interface are strongly reconstructed, with the number of 
holes in the Cu d 3z 2_ x states being close to the num- 
ber of Cu d x 2_ y 2 holes. A mechanism for such orbital 
reconstruction was proposed. 

In the present study, we focus our attention on the elec- 
tronic band structure of YBa2Cu307/Lai_ x Ca 2: Mn03 
superlattices and their XMCD and XAS spectra. We 
show that the magnetic circular dichroism at the Cu L2,3 
edges appears due to the hybridization between com- 
pletely filled Cu d 3z 2_ 1 states and unoccupied Mn d 3z 2_ x 
states split by the on-site exchange interaction. The num- 
ber of Cu d 3z 2_ 1 holes created by the hybridization is 
much smaller than the number of holes in the Cu d x 2_ y 2 
states. Nevertheless, even this tiny number of d 3z 2_ 1 
holes is sufficient to produce appreciable Cu L2,3 XMCD. 



The robustness of the results is verified by studying the 
influence of doping, atomic relaxation and correlation ef- 
fects on the electronic structure and XMCD spectra of 
YBCO/LCMO SL. 

The paper is organized as follows. The details of 
the calculations are described in Sec. |TTJ A structural 
model of YBCO/LCMO superlattices considered in this 
study is presented in Sec. Mil Section [TV] is devoted 
to the electronic structure and XMCD spectra of the 
YBCO/LCMO SL calculated using the local spin-density 
approximation (LSDA) and the LSDA+C7 method. Fi- 
nally, the results are summarized in Sec. [Vj 



II. COMPUTATIONAL DETAILS 

The electronic structure and XMCD spectra of 
YBCO/LCMO SL were investigated by means of ab 
initio calculations based on the spin-polarized density 
functional theory. The energy band structure, den- 
sities of states (DOS), and magnetic moments were 
calculated within the local spin-density approximation 
(LSDA) using the projected augmented plane wave 
(PAW) metho d 15 ' 16 ! 17 and the linear muffin-tin orbital 
(LMTO) method^ The cutoff energy of the plane-wave 
expansion was 500 eV and the spacing between k points 
was 0.03 A -1 . 

XMCD calculations were performed using the spin- 
polarized fully relativistic LMTO method (see Refs. [l9l l 20l 
for details) . To investigate the effect of electronic corre- 
lations on the XMCD spectra the rotationally invariant 
LSDA+J7 method 21,2 ™ with the double-counting term ap- 
proximated by the atomic limits was used. The screened 
on-site Coulomb repulsion U was considered as an ad- 
justable parameter. We used U=i to 5 eV for the Mn 
d states and C/ =7.5 eV for the Cu d states. For the 
exchange integral J the value of 1.0 eV estimated from 
constrained LSDA calculations was used. 

The finite lifetime of a core hole was accounted for 
by folding the XMCD spectra with a Lorentzian. The 
widths of the L2, and L 3 core level spectra, with 
T L2 =0.97 eV, r La =0.36 eV for the Mn and T L =1.04 eV, 
r 1,3=0. 61 eV for the Cu were taken from Ref. |24|. Finite 
instrumental resolution was accounted for by a Gaussian 
of 0.5 eV. 



III. STRUCTURAL MODEL 

The crystal structure of YBCO /LCMO SL was approx- 
imated by a tetragonal (P4/mmm) unit cell consisting 
of one YBCO and up to 3 LMO unit cells stacked along 
the c direction. The interface is formed by a BaO plane 
from the YBCO side and a Mn02 plane of LMO, which 
replaces a CuO plane of YBCO. The unit cell consisting 
of one YBCO and one LMO layer is shown in Fig. [T] 
For the in-plane lattice constants a=6=3.94 A the ex- 
perimental value for pseudo-cubic LMO is used. Thus, 
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FIG. 1: (Color online) Schematic representation of the crys- 
tal structure of the lxYBaCuO/lxLMO superlattice. The 
tetragonal unit cell used in the calculations is marked by the 
double-arrowed line. The number of electrons provided by 
each plane to Cu and Mn d shells is given in parenthesis. 



the orthorhombic distortion of bulk YBCO is neglected 
and the YBCO layer is slightly expanded in the ab plane. 
The c constant of 15.61 A is taken to be equal to the sum 
of the lattice constants of bulk YBCO (cybco=H-67 A) 
and pseudo-cubic LMO (clmo=3.94 A). 

Interfacial ions are subject to forces different from 
those in the bulk and, thus, can change their positions. 
The mere transfer of charge at an interface can lead 
to "electronic reconstruction". 25 Since the exact atomic 
structure at the YBCO/LCMO interface is not known, 
atomic positions were optimized using the PAW method. 
During the optimization the a and c lattice constants 
were fixed to the above mentioned values. As can be ex- 
pected, the most sensitive to the optimization are the po- 
sitions of the atoms which are close to the YBCO/LCMO 
interface. An oxygen ion denoted as Oea, which sits in 
the BaO plane between the nearest Cu and Mn ions, 
moves in the course of the optimization towards Mn. As 
a result, the Cu-Osa distance increases from 2.29 A in 
bulk YBCO to 2.37 A in the superlattice and, thus, be- 
comes even larger as compared to the separation of 1.98 
A between Cu and Oc u ions within the Cu02 plane. The 



Mn-Oea distance decreases from 1.97 A to 1.82 A. The 
optimized separation between Cu and Mn ions is slightly 
(4.18 A vs. 4.15 A) larger than between Cu ions resid- 
ing in Cu02 and CuO planes in bulk YBCO. The struc- 
ture optimization leads also to weak buckling of Mn02 
planes. Nevertheless, the Mii-Omh distance (1.98 A) re- 
mains practically the same as the Mn-OLa one (1.97 A). 
As the LaO plane is a mirror plane, the positions of OLa 
ions do not change during the optimization. Thus, the 
interfacial Mn ion is surrounded by a slightly distorted 
06 octahedron formed by Osa, OLa, and 4xOmii ions. 

Atomic positions were also optimized for 
lxYBC0/2xLM0 SL with c=19.55 A (c Y bco + 2c LMO ) 
containing 2 LMO unit cells. In this case, the Cu-Osa 
and Mn-Oea distances of 2.48 A and 1.80 A were 
obtained, which gives the separation of 4.28 A between 
Cu and Mn ions at the interface. As the c lattice 
constant was fixed, the expansion at the interface occurs 
at the expense of some decrease of the distances between 
Cu02 planes and between the "bulk" and "interfacial" 
Mn02 planes. 

Under the assumption that O p states are completely 
filled and s and d states of Ba, Y, and La are empty, Mn 
and Cu d states in the SL shown in Fig. Q] are filled with 
24 electrons. Taking into account that the unit cell con- 
tains 2 Cu and 2 Mn ions and assuming that some charge 
5 may be transferred across the interface, the occupations 
of the d shells are Cu d 9±s and Mn d 3zfS . Hence, the oc- 
cupation of Mn d states in this model SL corresponds 
to Lai-^Ca^MnOa with x ~ 1. With each LMO layer, 
consisting of one Mn02 and one LaO plane, added to a 
model unit cell, the number of electrons which fill the d 
states of transition metal ions increases by 4, while the 
number of Mn ions only by 1 . Consequently, the average 
occupation of Mn d states increases with the increase of 
the LMO layer thickness. This should be taken into ac- 
count when comparing the results calculated for different 
model YBCO/LMO SL. 



IV. RESULTS AND DISCUSSION 
A. XMCD spectra 

The x-ray magnetic circular dichroism in the 
YBCO/LCMO SL has been measured at the Mn and 
Cu L2,3 edges by Chakhalian et al. (Ref . [Io|) . The mea- 
surements were carried out at 30 K in an applied field 
of 500 Oe on a superlattice of alternating 100 A-thick, o 
axis-oriented YBa2Cu307 and La2/3Cai/3Mn03 layers. 
A circularly polarized x-ray beam was shined at a small 
angle to the film surface, so that the x-ray propagation 
direction was almost parallel to the LCMO magnetiza- 
tion oriented in the ab plane. 

Figure [2] shows XMCD spectra at the Mn and Cu L2,3 
edges calculated within LSDA for the lxYBCO/lxLMO 
SL together with the experimental spectra from Ref. [Tol . 
The corresponding normalized theoretical x-ray absorp- 
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FIG. 2: (Color online) The theoretical and experimental 
(Ref. 03) XMCD spectra of the YBCO/LCMO superlattice 
at the Mn and Cu Z/2,3 edges. The experimental Cu Z/2,3 
XMCD spectrum is multiplied by a factor of 2. The insets 
show theoretical XAS spectra averaged over two circular po- 
larizations (solid line). The Cu 1/2,3 XAS calculated with the 
x-ray polarization along the c axis is plotted by dashed line. 
All the spectra are plotted relative to the corresponding L3 
absorption edge marked by vertical dashed lines. 



tion spectra averaged over two circular polarizations are 
plotted by solid lines in the insets in Fig. [2l The LSDA 
calculations reproduce very well the experimental Cu 
Z/2,3 XMCD spectra, although the calculated magnetic 
dichroism has a small oscillatory structure at the high 
energy part of the Cu L2.3 spectra which is not detected 
by the experiment probably due to a weak signal. XMCD 
peaks at the Mn and Cu L3 edges are of opposite signs, 
as it is expected when the magnetization induced on Cu 
is antiparallel to Mm d . The theory overestimates the 
dichroism at the Cu 1/2,3 edges approximately by a fac- 
tor of two in comparison with the experimental observa- 
tion. This can be explained by the fact that our calcula- 
tions were performed for the unit cell containing a single 
YBCO cell, with the Cu ions in both Cu02 planes polar- 
ized due to their proximity to magnetic LCMO layers. In 
the experimental samples, on the other hand, YBCO lay- 
ers, consisting of about 3 unit cells of bulk YBa2Cu307, 
are significantly thicker. It is natural to expect that 
only Cu ions in the interfacial Cu02 planes have a non- 
zero component of the magnetic moment along the light 
propagation direction and, consequently, contribute to 



the measured XMCD signal, whereas inner Cu02 planes 
are XMCD-inert. With the increase of the YBCO layer 
thickness in the model unit cell, the XAS intensity would 
increase proportionally to the number of CuO planes, but 
not the dichroism. Since the calculated XMCD spectra 
were normalized by the intensity of the corresponding 
x-ray absorption spectra, the normalized XMCD signal 
would tend to decrease with the increase of the YBCO 
layer thickness, which would improve the agreement be- 
tween the calculations and the experiment. 

Our calculations for the lxYBCO/lxLMO SL result 
in strong linear dichroism at the Cu £2,3 edges. The XAS 
spectrum calculated for the x-ray polarization e parallel 
to the c axis is compared to the spectrum averaged over 
two circular polarizations, which in this geometry corre- 
sponds to the averaging over e\\b and e\\c polarizations, in 
the inset to the lower panel of Fig. [2l A prominent peak 
at ^1 eV in the averaged spectrum is strongly suppressed 
when e becomes parallel to c. The strong XLD at the Cu 
£2,3 edges is at variance with the interface-sensitive TEY 
measurements^ which showed a very weak polarization 
dependence of Cu L2.3 XAS. 

In the case of Mn, the experimental L2.3 XMCD spec- 
trum shows two additional positive peaks at the energies 
around 4 eV and 10 eV which are absent in the theoretical 
spectrum. The cause of such a disagreement is yet un- 
clear. One of possible reasons may be many-body effects, 
e.g., the interaction of the partially filled Mn d states with 
a 2p core hole, which is not taken into account in the the- 
oretical calculations. The Mn L2.3 XMCD spectrum may 
change if the magnetization in interfacial Mn02 planes 
is reversed relative to the rest of the LCMO layer or it 
may be strongly influenced by possible orbital ordering 
in LCMO. These questions need further theoretical con- 
sideration. 

In the following we focus on the origin of the mag- 
netic circular dichroism at the Cu L2.3 edges, but first 
we turn our attention to the electronic band structure of 
YBCO/LMO superlattices. 



B. Energy band structure 

Symmetry-resolved densities of Cu d, Mn d, and 
O p states calculated within LSDA for the model 
lxYBCO/lxLMO SL (Fig. [TJ assuming ferromagnetic 
ordering of Mn magnetic moments are presented in Fig. [3] 
Mn d states are split by the octahedral component of 
the ligand field into rather narrow d xy and d xz ^ yz ("£23") 
states and much wider d x 2_ y 2 and c?3 Z 2_ 1 ("e s ") ones. 
The latter are strongly hybridized with, respectively, 
O m n Px,y and Ojja and OLa Pz states. The d^ z 2_i states, 
or more precisely their antibonding combinations with O 
p states, are shifted to higher energies due to the short- 
ening of the Mn-Oea bond at the interface. Because of 
the small thickness of a LMO layer, they are narrower 
than the d x 2_ y 2 states. 

The on-site Hund's coupling splits the Mn t2 g states 
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FIG. 3: (Color online) Spin- and symmetry-resolved LSDA 
densities of states, N(E), for lxYBCO/lxLMO SL. The 
Fermi energy is at zero. 



into almost completely occupied majority- and unoccu- 
pied minority-spin states. The t2 g states provide the 
dominant contribution (2.21/ib) to the magnetization 
-^Mn=2.75/XR inside a Mn sphere. The rest comes from 
the Mn d x 2_ y 2 (0.32/ib) and d 3z 2_ 1 (0.21^B) states which 
are also split by the on-site exchange interaction. While 
the minority-spin e g states are empty, the bottom of the 
majority-spin d^ z 2_i and, especially, d x i_ y i states is par- 
tially filled. 

Figure [3] shows that, as it is usual for cuprates, all Cu 
d states, except for d x i_ y i ones, lie well below the Fermi 
level (sf)- The d x 2_ y 2 states, interacting strongly in the 
ab plane due to dpa-type bonds with Oc u Px,y states, 
form 3 eV-wide bands in the range from —1 to 2 eV, as 
can be seen in Fig. 0] which shows the contributions of the 
<^3z 2 -i an d d x 2_ y 2 states of Cu and Mn to the bands cal- 
culated for the lxYBCO/lxLMO SL. The bands origi- 
nating from Cu and Mn d x 2_ y 2 states have similar disper- 
sions. The Cu d x 2_ y 2 bands for both spin directions are 
nearly degenerate and cross Ef- The majority-spin Mn 
bands also cross e f , whereas the minority-spin ones lie 2 
eV higher and are completely empty. Mn <i 3z 2_ 1 -derived 
bands are split into bonding and antibonding ones by 
strong interaction between the d 3z 2_i orbitals of Mn ions 
from two Mn02 planes via OLa Pz states. The bonding 
majority-spin Mn <i 3z 2_ 1 bands are found just above ef- 

Finally, Cu d 3z 2_ 1 bands for both spin directions lie 
about 2 eV below Ef and hybridize strongly with O p 
bands. Owing to the hybridization of the Cu d^ z 2_^ 
states with the spin-split Mn d 3z 2_ 1 states, which is me- 
diated by Oea Pz orbitals, the densities of the majority- 
and minority-spin Cu d^ 2 -!, as well as Oe a Pz, states 



are remarkably different (see Fig. [3]), with the differences 
being mostly restricted to the occupied bands. However, 
due to the Cu <i3 Z 2_i-OBa p z -Mn d 3z 2_ 1 hybridization, 
the Cu d^ z 2_i states contribute also to the unoccupied 
Mn d 3z 2_i -derived bands. Because of the large (2.37 A) 
Cu-Oea separation and a large energy difference between 
the Cu and Mn d 3z 2_ 1 states, this contribution is quite 
small and can hardly be seen in Fig. [3] It is, however, im- 
portant for understanding the origin of Cu L2.3 XMCD 
spectra presented in Fig. 

It is well documented that LSDA fails to produce a 
local Cu magnetic moment in HTSC cupratesi 23 i 26 In 
the model lxYBCO/lxLMO SL, a tiny Cu moment 
Mcu=— 0.002/xb, i.e., directed oppositely to the Mn mo- 
ment, is induced by the proximity of a ferromagnetic 
Mn02 plane. The dominant contribution of — 0.014/xb 
to Mcu is provided by the d3 Z 2_! states. It is, however, 
partially compensated by the magnetization of opposite 
sign (0.012/is) induced in the d xz . yz chanel. The con- 
tributions of the d xy (0.003/ib) and d x 2_ y 2 (— 0.003/ib) 
states are even smaller and cancel each other. Such a 
small difference in the occupations of the majority- and 
minority-spin Cu d states produces negligibly small ex- 
change splitting of the d shell. The d% z 2_i states become 
spin-polarized as a result of their hybridization with the 
spin-split unoccupied Mn d 3z 2_ 1 states. The majority- 
spin Mn d^ z 2_ 1 bands disperse just above £p, about 3 
eV higher than the center of Cu d$ z 2_ 1 states, while the 
minority-spin bands are shifted by the on-site exchange 
splitting A ex ~ 2 eV to higher energies (Fig. 2]). Be- 
cause of this energy difference, the Cu <i3 Z 2_i-OBa Pz~ 
Mn da z 2 -i hybridization creates a larger number of holes 
in the majority- than in the minority-spin Oea Pz and Cu 
^3z 2 -i states, as illustrated schematically in Fig.® Con- 
sequently, the sign of the corresponding contributions to 
Cu and Oea magnetizations is opposite to the magne- 
tization of Mn. The negative Oe a Pz spin-polarization 
is, however, over-compensated by the positive one, in- 
duced by the hybridization of p XtV orbitals with Mn d xz ,yz 
states: as the majority-spin d XZtVZ states are completely 
filled, only minority-spin p XtV holes are created and the 
calculated total Os a spin moment of 0.13/zb is parallel to 
M M „. 

The spacial distribution of the spin-density (SD) calcu- 
lated for the lxYBCO/lxLMO SL is presented in Fig.® 
Large positive iso-surfaces surrounding Mn ions (right 
panel of Fig. [6]) depict the spin-density of the completely 
spin-polarized Mn t2 g states. The Mn <i £CZj j /z -OBa Px,y 
hybridization is responsible for the positive SD surface 
around Ob&- The left panel of Fig. [5] shows clearly that 
negative SD surfaces surrounding Oea and Cu ions ap- 
pear due to the polarization of Oua Pz and Cu d 3z 2_ 1 
states, respectively, which, as it is explained above, is 
caused by their interaction with the spin-split unoccu- 
pied Mn d3 Z 2_! states. 

Test calculations in which Oea was shifted along the 
line connecting Mn and Cu ions, while all other atomic 
positions were kept fixed, showed that the calculated Cu 
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FIG. 4: (Color online) The energy band dispersion in the fat-band representation calculated within LSDA for the model 
lxYBCO/lxLMO SL shown in Fig. [1] The size of a circle is proportional to the weight of the corresponding state in a Bloch 
wave function at a given k-point. 
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FIG. 5: (Color online) A sketch of the hybridization between 
the occupied Cu d 3z 2_ 1 states and the unoccupied Mn d 3z 2_ 1 
states split by the on-site exchange interaction A ex . The 
contribution of the Cu d 3z 2_ 1 orbitals to bonding (B) and 
antibonding (AB) states is shown by thick lines. The nega- 
tive spin-polarization of the Cu d 3z 2_ 1 states appears due to 
their larger weight in the majority-spin AB states than in the 
minority-spin ones. 

moment is very sensitive to the position of the interfa- 
cial Oea ion. The variation of the Mn-Oea and Cu- 
Orja distances has a twofold effect: first, the increase 
of the Mn-Osa separation weakens the Mn d 3z 2_ 1 -Osa 
p z hybridization, which leads to a downward shift of 



the mostly unoccupied d 3z 2_i bands relative to d x 2_ y 2 
ones. The majority-spin d 3z 2_ 1 states become more pop- 
ulated, which provides the dominant contribution to the 
increase of the Mn spin moment from 2.75/zb, when the 
optimized Oe a position is used (dMn-o Ba = l-82 A), to 
3.01/ib, when C-Ba is equidistant to the Mn and Cu sites 
(dMn-o Ba = 2-09 A). It should be noted, that in spite 
of the increase of the Mn-Oea distance, the negative 
spin-polarization of the p z states increases as Oea moves 
away from Mn. This can be explained by the decrease 
of the energy separation between the majority-spin Mn 
d^ z 2_i and Oea Pz states as a result of the downward 
shift of the former and an upward energy shift of the 
latter. Second, strengthening of the Cu d 3z 2_ 1 -Osa Pz 
hybridization due to the decrease of the Cu-Osa distance 
increases the Cu c?3 Z 2_ 1 weight in the Mn d 3z 2_ 1 -derived 
bands and, consequently, the spin-polarization of the Cu 
^3z 2 -i states. As Os a is shifted from the optimized po- 
sition (dcu-o B a— 2.37 A) to the middle of the Cu-Mn 
line (dcu-o Ba =2-09 A), Mq u increases from —0.002 to 
-0.036/LtB- 



C. Origin of XMCD at Cu L 2 ,3 edges 

The calculated Cu L2 and L3 XMCD spectra and the 
spin-polarization of Cu d 3z 2_ 1 and d x 2_ y 2 states, i.e., the 
difference of the densities of corresponding majority- and 
minority-spin states, are plotted together in lower and 
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FIG. 6: (Color online) The calculated spacial distribution of 
the spin-density in the lxYBCO/lxLMO super-lattice. Left 
and right panels show negative and positive const ant- spin- 
density iso-surfaces, respectively. The cut is through a [100] 
plane, with z axis parallel to c. 



upper panels of Fig. [7J In order to make the similarities 
more clear, the spectra are plotted without the life-time 
and instrumental broadening and arc aligned to the cor- 
responding absorption edges which are shifted to zero 
energy. Figure [7J shows unambiguously that the energy 
dependence of the calculated XMCD signal at the Cu 
and L3 edges is determined by the spin-polarization of 
the Cu d 3z 2_ 1 states. 

Qualitative understanding of the Cu L2.3 XMCD ab- 
sorption spectra is provided by the analysis of the matrix 
elements for 2p 1 / 2 ,3/2 - * 3e? 3 / 2 ,5/2 transitions. XMCD at 
the L2 } 3 edges is mostly determined by the strength of 
the spin-orbit coupling of the initial 2p core states and 
spin-polarization of the final 3d states while the exchange 
splitting of the former as well as the SO coupling of the 
latter are usually of minor importance. 27 Squared angu- 
lar matrix elements for the absorption of a photon with 
left (p=+l) or right (p—— 1) circular polarization, fij_, 
by exciting an electron from the SO-split Pi/2,3/2 core 
levels to spin- and crystal-field split d states are collected 
in Tables Q] and [TTJ The matrix elements were calculated 
in the dipole approximation neglecting the SO splitting 
of the d states and assuming that x-rays propagate along 
the spin quantization axis which is parallel to z. 
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FIG. 7: (Color online) The comparison of the calculated Cu 
1/2 and L3 spectra (lower panel) and the difference of the 
densities of the majority- and minority-spin Cu d 3z 2_ t and 
d x 2_ y 2 states (upper panel). The unbroadened spectra are 
plotted relative to the corresponding absorption edges. 



TABLE I: Squared angular matrix elements Q,± for dipole 
allowed transitions at the L2 edge from an initial core state 
with the total angular momentum j = 1/2 and its projection 
rrij to a final 3d valence state (|/)) with the absorption of 
left (p = +1) or right (p = —1) circularly polarized x-rays. 
ft± are multiplied by 90. The intensity of unpolarized x-ray 
absorption and XMCD spectra are proportional to f2+ + fl 2 ^ 
and f2+ — ftL , respectively. 





P = 


+ 1 


p = 


-1 






1/) 


m,j 


n 2 + 


mj 


nl 




n 2 + - n 2 _ 


dfyz 


-1/2 


3 


-1/2 


3 


6 





dfzx 


-1/2 


3 


-1/2 


3 


6 





d^xy 


1/2 


12 






12 


12 








1/2 


4 


4 


-4 


d-\ x 2 — y 2 


1/2 


12 






12 


12 


dlyz 


1/2 


3 


1/2 


3 


6 





dlzx 


1/2 


3 


1/2 


3 


6 





dixy 






-1/2 


12 


12 


-12 


dl3z 2 -l 


-1/2 


4 






4 


4 


d lx 2_ y 2 






-1/2 


12 


12 


-12 



8 



TABLE II: Squared angular matrix elements flj. for dipole 
allowed transitions at the Lz edge. The notations are the 
same as in Table [I] 
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By comparing Tables Q] and [TT1 we can conclude that 
neglecting the radial matrix elements and different life- 
times of 2p\/2 and 2p 3 / 2 core holes, the unpolarized XAS 
intensity, proportional to Vt 2 ^ + fi 2 ,, is twice as high at 
the L3 edge as at the L2 one. The corresponding XMCD 
signals, Vl 2 + — fi 2 ,, are, however, of the same magnitude, 
with the sign of the dichroism at the L 2 edge being oppo- 
site to the Lz one. This can be clearly seen in the lower 
panel of Fig. [7] where the unbroadened Cu L2 and L3 
XMCD spectra are plotted. The almost twofold differ- 
ence in the intensities of the L2 and L3 XMCD spectra 
presented in Fig. [5] appears because of the large differ- 
ence in the inverse life times of 2pi/2 (Fl 2 =1.04 eV) and 
2p 3 /2 (r L3 =0.61 eV) core holes,2£ 

The dependence of the matrix elements on the symme- 
try of the hnal 3d state can be summarized as follows: the 
dominant contributions to the x-ray absorption and to 
the dichroism, i.e., the difference of the absorption prob- 
abilities for left and right polarized photons, fl 2 ^ — f2?_, 
come from the transitions to d xy and d x 2_ y 2 states. The 
transitions to d 3z 2_i states are weaker and their contri- 
bution to the dichroism is of the opposite sign. Finally, 
the probabilities to excite a core 2p electron to a d zx or 
d yz state with either left or right polarized photons are 
equal so that these transitions do not contribute to the 
XMCD spectra. 

It should be emphasized that the contributions to the 
dichroism from transitions to final states with differ- 
ent spin projections have opposite signs. Consequently, 
if majority- and minority-spin d states with particular 
orbital character are degenerate their contributions to 
XMCD spectra cancel each other. This is exactly what 
happens to the Cu d x 2_ y 2 contribution to the XMCD 
spectra calculated for the lxYBCO/lxLMO SL: al- 
though the density of the unoccupied Cu d x 2_ y 2 states 
is much higher than of the d 3z 2_ 1 ones (see Fig. [3]) and, 
consequently, transitions to the final states of the d x 2_ y 2 
character control the x-ray absorption at the L2 and L3 



edges, the majority- and minority-spin d x 2_ y 2 bands are 
degenerate and their contribution to the dichroism is very 
weak. 

In contrast to d x 2_ y 2 , the unoccupied Cu d 3z 2_i states, 
which, as explained in Sec. IIVB1 follow the dispersion of 
the spin-split Mn d 3z 2_ 1 bands, are strongly polarized 
due to the Cu c^^^-Osa p z -Mn rf 3z 2_ 1 hybridization. 
As the minority-spin Mn bands are shifted to the energies 
above 2 eV by the strong exchange interaction within the 
Mn d shell (see Fig. 2]), only the majority spin Cu d^ 3z 2_ 1 
holes are created in the energy range up to 2 eV above 
£p. Although the number of the df 3z 2_ 1 holes is very 
small, it is sufficient to produce the appreciable XMCD 
signal at Cu L 2: 3 edges. 

According to Tables U and |TT] the matrix elements for 
the transitions to the final states of the d 3z 2_ 1 charac- 
ter arc small. It should be noted, however, that they 
were derived for the light propagation and magnetiza- 
tion directions parallel to z, i.e., to the tetragonal c axis, 
whereas the XMCD spectra shown in Fig. [5] were calcu- 
lated for the experimental geometry^ in which circularly 
polarized x-rays propagate along the direction of the Mn 
magnetization which is parallel to the a axis. The de- 
pendence of the XMCD spectra on the light propagation 
direction can be estimated by expanding the Cu d% z 2_i 
orbital in terms of cubic harmonics defined in the rotated 
frame with z'\\a: 

A l A' ^^A< 

dz? = ~ 2 22 + — d x 2 -y 2 ■ 
Because of the appearance of the d' 2 , 2 contribution in 

x y 

the expansion of the corresponding wave functions, the 
Cu i 2 ,3 XMCD spectra calculated in the z'\\a geometry 
should change sign and become more intense compared 
to the z\\c spectra. This strong orientational dependence 
of the Cu XMCD spectra was verified by performing cal- 
culations for the lxYBCO/lxLMO SL with M Mn ||e. In- 
deed, due to the weakness of the matrix elements for the 
transitions to the d 3z 2_ 1 -like final states, the magnitude 
of the Cu L2,3 XMCD spectra calculated in this geome- 
try is smaller and their signs are opposite compared to 
the spectra shown in Fig. [5] which were obtained with the 
light propagation and Mu n directions parallel to a. 



D. Robustness of calculated Cu L2,3 XMCD 

The proposed explanation for the strong XMCD at the 
Cu £2,3 edges and robustness of the calculated spectra 
have been confirmed by a number of additional calcula- 
tions. As explained in Sec. IIVBI the Cu magnetization 
induced by the proximity to the LMO layer is sensitive 
to the position of Oea- This is also true for the XMCD 
spectra. The increase of the Cu <i3 Z 2_i-OBa Pz hybridiza- 
tion strength with the decrease of the Cu-Oea distance 
leads to the growth of the Cu d3 Z 2_i weight in the Mn 
d| 3z 2_ 1 -derived bands. Because of the increased number 
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of Cu d|3 Z 2_i holes the calculated Cu L2.3 XMCD spec- 
tra become more intense as Oea shifts from Mn toward 
Cu. 

As was mentioned above the number of electrons oc- 
cupying Cu (d 9±l5 ) and Mn (d 3 ^) 3d shells in the model 
lx YBCO/lxLMO SL corresponds to the composition 
of the Lai-^CazMnOa layer with x ~ 1. By using the 
virtual crystal approximation we performed also LSDA 
calculations with the increased number of valence elec- 
trons corresponding to Ca doping of x=2/3 and 1/3. The 
added electrons fill mainly Cu d x 2_ y 2 and majority-spin 
Mn d x 2_ y 2 bands crossing the Fermi level. Due to the in- 
crease of the Fermi energy, the Mn d| 3z 2_ 1 -derived bands 
moves closer to Ep but their occupation does not change. 
As a consequence, the change of the electron count has 
only a very weak effect on the Cu XMCD spectra. 

Finally, the effect of the increase of the LMO layer 
thickness on the Cu L 2 ,3 XMCD spectra was studied 
by performing LSDA calculations for lxYBCO/nxLMO 
(rt=2, 3) supcrlattices with FM alignment of Mn mag- 
netic moments. A LMO layer in these SL consists 
of 2 and 3 unit cells of pseudo-cubic LaMn03, re- 
spectively, i.e., 3xMn0 2 +2xLaO and 4xMn0 2 +3xLaO 
planes. The distance between the Mn0 2 planes was 
fixed to clmo— 3-94 A. For ions at the interface the po- 
sitions obtained from the structure optimization for the 
lx YBCO/lxLMO SL were used. In contrast to the lat- 
ter, in which both M11O2 planes are at the interface with 
the YBCO layer, the thicker superlattices have addition- 
ally one (n=2) or two (n=3) "bulk" Mn0 2 plane(s). Mn 
ions in these bulk Mn0 2 planes are surrounded by an 
undistorted Oq octahedron, with all the Mn-0 distances 
of 1.97 A being equal to that in pseudo-cubic LaMn03. 

Assuming that the Cu valency does not depend on the 
LMO layer thickness, the effective Mn d electron count 
in these SL increases compared to the SL with n=l and 
corresponds to Ca content of £=2/3 for n=2 and x=l/2 
for n—i. In order to separate the effects of the LMO layer 
thickness and that of the varying occupation of the Mn 3d 
shell, XMCD calculations for all three SL were performed 
for the same effective Mn d occupation corresponding to 
x=2/3. 

The magnetic moment of ~3.2/iB calculated for the 
bulk Mn ions is higher than for the interfacial ones 
and close to the LSDA value obtained for pseudo-cubic 
Lai_2;Ca x Mn03 with the same 2=2/3. The moment of 
the interfacial Mn ions and the magnetization induced on 
Cu remain nearly constant independently of n. As a con- 
sequence, the Cu L 2i3 XMCD spectra calculated for the 
SL with ro=l, 2, and 3 are nearly identical, which shows 
that they are not sensitive to the LMO layer thickness. 



E. Effect of electronic correlations 

The influence of electronic correlations on the Cu L 2j 3 
XMCD spectra calculated for the lx YBCO/lxLMO 
SL was investigated using the rotationally invariant 



LSDA+J7 method. 21,22 In order to understand the 
LSDA+f/ results it is worth recalling that the expres- 
sion for the additional orbital dependent LSDA+/7 po- 
tential Vi acting on i-th localized orbital becomes particu- 
larly simple if non-spherical contributions to the screened 
Coulomb U and exchange J integrals are neglected: 

Vi = QJ-J)(\-m), 

where m is the occupation of the i-th orbital. Thus, 
the occupied (unoccupied) states are shifted downwards 
(upwards) by ~ (U — J)/2. 

When U is applied only to the Mn d states its main 
effect is to increase the splitting between the occupied 
majority- (n,; ~ 0.9) and unoccupied minority-spin (n, ~ 
0.1) Mn "t 2g " states. This causes some increase of Mmh 
from the LSDA value of 2.8/zb to 3.0/ib for U—3 eV and 
3.3/iB for U—5 eV. The minority-spin d 3z 2_ 1 and d x 2_ y 2 
states with n, ~ 0.2 are also shifted to higher energies. 
The occupation numbers of the majority-spin Mn d 3z 2_ 1 
and d x 2_ y 2 states are close to 0.5 because of their strong 
hybridization with O p states and partial filling of the 
corresponding bands (see Fig. S}. As a consequence the 
majority-spin Mn d 3z 2_ 1 - and d x 2_ y 2 -derived bands are 
almost unaffected by U. Hence, the Cu Z/2,3 XMCD spec- 
tra, which are mostly determined by the Cu d^z 2 -i~ Mn 
d| 3z 2_x hybridization, show weak dependence on U ap- 
plied to the Mn d states. 

LSDA+C7 calculations with U—7.5 eV acting within 
the Cu 3d shell were performed for an \[2a x \/2a x c 
unit cell of the lxYBCO/lxLMO SL. In this enlarged 
cell two Cu sites in each Cu0 2 plane become inequiv- 
alent which enables calculations with antifcrromagnetic 
(AFM) ordering of Cu moments. The calculations were 
performed for a collinear magnetic arrangement with 
both Cu and Mn magnetic moments parallel to the 
tetragonal a axis. 

LSDA calculations for the double cell converge to the 
same magnetic solution which was obtained for the single 
cell and is characterized by FM ordering of Mn moments 
and a tiny FM moment induced across the interface on 
Cu ions. However, when [7=7.5 eV is applied to the 
Cu d states and AFM ordering within a Cu0 2 plane is 
allowed, Cu ions acquire a spin moment of ~0.4/xb- The 
moment appears due to splitting of the partially occupied 
majority- and minority-spin Cu d x 2_ y 2 states by the on- 
site Coulomb repulsion. The other Cu d states, including 
^3z 2 -i ones, which are completely filled already in LSDA 
calculations, are shifted ~3 eV downwards to the bottom 
of O p bands. 

Because of the proximity to FM Mn0 2 planes the Cu 
moments which are parallel (0.44/ib) and antiparallcl 
(— 0.45//b) to Myin are slightly different, with the dif- 
ference being comparable to the LSDA moment induced 
in the Cu d 3z 2_ 1 chanel by the Cu d 3z 2_ 1 -OB a p 2 -Mn 
d 3z 2_ 1 hybridization as discussed in Sec. UVBl Although 
the hybridization causes a small (^O.OI^b) difference in 
the moments of two inequivalent Mn ions which have Cu 
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FIG. 8: (Color online) Cu L 2 ,3 XMCD spectra calculated 
using the LSDA+C7 method with [7=7.5 eV applied to the 
Cu 3d states for Cu ions with the magnetic moment parallel 
(dashed line) and antiparallel (dash-dotted line) to Mytn- The 
averaged spectrum is plotted by solid line. For comparison, 
the LSDA spectrum from Fig. [2] is also shown (dotted line). 
The spectra are plotted relative to the L3 absorption edge. 

neighbors with either parallel or antiparallel magnetiza- 
tion, the calculated Mm d of 2.75z^b and 2.76zxb are very 
close to the LSDA value. 

Figure [8] shows that the spin-polarization of the Cu 
d x 2_ y 2 states affects strongly the £2,3 XMCD spectra 
calculated for Cu ions with the moment parallel (Cuj) 
and antiparallel (CujJ to Mmh- As the majority- and 
minority-spin Cu d x 2_ y 2 states are split by the on-site 
Coulomb repulsion, transitions to the d x 2_ y 2 final states 
give a substantial contribution to the dichroism. In the 
geometry with the x-ray propagation and magnetization 
directions parallel to the c axis the contribution of the 
d x 2_ y 2 states would be the dominant one and the Cuj 
and Cu| XMCD spectra would be of the same shape but 
opposite signs. Since the calculations were performed 
with Mcu||a the matrix elements for the transitions to 
the d x 2_ y 2 states are strongly reduced and their contri- 
bution is comparable to the d 3z 2_ 1 one, which is induced 
by the hybridization with FM ordered Mn and, thus, 
is of the same sign for both Cu ions. As a result, the 
shapes of the Cu-p and Cuj L3 XMCD spectra are dif- 
ferent (Fig. |5J), whereas the Li spectra are of the same 
sign although the Cu-f and Cuj moments are antiparal- 
lel. Nevertheless, when the individual Cu-f and Cux spec- 
tra are averaged, the net Cu XMCD spectrum becomes 
similar to the LSDA spectrum calculated for the single 
lxYBCO/lxLMO cell, which is also shown in Fig. for 
comparison sake. The magnitude of the averaged spec- 
trum is two times smaller, which may be explained by 
decrease in the weight of the Cu d 3z 2_ x states in the 
minority-spin Mn d 3z 2_i bands caused by their down- 
ward shift in the LSDA+[7 calculations. Hence, even 
in the presence of the sizable Cu moment the dichroic 
signal at the Cu L 2 ,3 edges is still governed by the spin- 



polarization of the Cu d 3z 2_! states induced by the Cu 
d 3z 2_ 1 -Mn d 3z 2_x hybridization via the p states of inter- 
facial Oe a ions. 

Note, that if the AFM aligned Cu d x 2_ y 2 moments 
were oriented perpendicular to the Mn magnetization 
and, consequently, to the x-ray propagation direction 
they would not contribute to the dichroism, which also in 
this case would be determined by the transitions to the 
Cu d 3z 2_i states. 



V. CONCLUSIONS 

We studied the electronic structure, XAS and XMCD 
spectra of YBCO/LMO superlattices by means of the ab 
initio PAW and fully-relativistic spin-polarized LMTO 
methods. 

Our results contradict to the orbital reconstruction 
model, proposed in Ref . 1 1 in order to explain weak XLD 
at the Cu £2,3 edges measured in the TEY mode. In 
all the calculations performed in the present work the 
number of Cu d 3z 2_i holes is significantly smaller than 
the number of d x 2_ y 2 ones. The large difference in the 
number of d$ z 2_i and d x 2_ y 2 holes causes strong linear 
dichroism of the calculated Cu L 2 ,3 XAS. 

On the other hand, the LSDA calculations reproduce 
very well the experimental Cu L 2 ,3 XMCD spectra of 
Ref. 0. We have shown that the theoretical XMCD 
spectra are proportional to the difference of the densi- 
ties of unoccupied majority- and minority-spin Cu d% z 2_i 
states. Although the Cu d 3z 2_ 1 -derived bands lie ~2 eV 
below the Fermi level, a small number of d$ z *—\ holes is 
created because of the Cu d^a^-OBa p z -Mn d 3z 2_i hy- 
bridization, which induces a small Cu magnetic moment 
opposite to the moment of Mn ions. However, this tiny 
number of d 3z 2_ 1 holes is sufficient to produce substantial 
XMCD at the Cu L 2 ,3 edges. 

The calculated Cu £2,3 XMCD spectra are very ro- 
bust: they are not sensitive to Ca doping of the 
Lai_ 3; Ca a ;Mn03 layer simulated using the virtual crystal 
approximation, nor to the LCMO layer thickness. Re- 
markably, the shape of the XMCD spectra obtained from 
LSDA+C/ calculations does not change much even when 
U=7.5 eV is applied to the Cu d states and AFM order 
with the Cu moment of 0.4/ljb is induced in the Cu02 
plane. 
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